SPIKES AND DIFFUSION WAVES 
IN ONE-DIMENSIONAL MODEL OF CHEMOTAXIS 
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Abstract. We consider the one-dimensional initial value problem for the viscous trans- 
port equation with nonlocal velocity u t = u xx — (u(K' * u)) with a given kernel K' £ 
We show the existence of global-in-time nonnegative solutions and we study 
their large time asymptotics. Depending on K 1 , we obtain either linear diffusion waves 
(i.e. the fundamental solution of the heat equation) or nonlinear diffusion waves (the 
fundamental solution of the viscous Burgers equation) in asymptotic expansions of solu- 
tions as t — > oo. Moreover, for certain aggregation kernels, we show a concentration of 
solution on an initial time interval, which resemble a phenomenon of the spike creation, 
' typical in chemotaxis models. 

m 

<^ ! 1. Introduction 

In this work, we study the large time behavior of solutions to the one- dimensional initial 
value problem 

(1.1) u t = u xx - (u(K' * u)) x for iGl, t> 0, 

(1.2) u{x, 0) = u (x) for x G R, 

^ where the aggregation kernel K' G L X (R) is a given function (the symbol "*" denotes 

the convolution with respect to the variable x) and the initial datum it G L 1 (R) is 
nonnegative. Such models have been used to describe a collective motion and aggregation 
phenomena in biology and mechanics of continuous media. In this case, the unknown 
function u = u(x,t) > is either the population density of a species or the density of 
particles in a granular media. The kernel K' in (II. ip can be understood as the derivative 
of a certain function K, that is, K' stands for dK/alx. We use this notation to emphasize 
that cell interaction described by equation (II. ip takes place by means of a potential K 
^ ■ (see [7j for derivation of such equations). 

Equation (II .ip and its generalizations, considered either in the whole space or in a 
bounded domain, have been studied in several recent works. First, notice that in the 
particular case of K(x) = e~' x '/2, equation (11.11) corresponds to the one-dimensional 
parabolic-elliptic system of chemotaxis 

(1.3) u t = Uxx - {uv x ) x , -v xx + v = u, x G R, t > 0. 
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Indeed, since K(x) = e~' x '/2 is the fundamental solution of the operator —d 2 x + Id, one 
can rewrite the second equation in (11.31) as v = K * u. Now, it suffices to substitute this 
formula into the first equation in (11.31) to obtain (11.11) . We refer the reader to the works 
[21 IH Q2J [6j [T71 [22] (this list is far from being complete) for mathematical results and for 
additional references on systems modeling chemotaxis. 

Next, one should mention the inviscid aggregation equation u t + V • {uVK * u) = 0, 
describing the evolution of a cell density, which was derived as a macroscopic equation 
from the so-called "individual cell-based model" [29, 7\. Here, the reader is referred to 
P, El El HE] for recent results on the existence and the blowup of solutions to the initial 
value problem for the inviscid aggregation equation. 

To handle diffusion phenomena, equations describing aggregation are supplemented 
with additional terms. One possible approach is to add a nonlinear term modeling a 
degenerate diffusion as in the porous medium equation, see e.g. [32J. Results and other 
references on the chemotaxis model with a degenerate diffusion can be found in [301 [5] 
and on more general aggregation equations in [T51 EU I2"U] . 

In several cases, the mechanism of spreading out of organisms resembles a Levy flight, 
hence, the anomalous diffusion is better modeled by nonlocal pseudodifferential operators. 
Recent works [21 HI [HI [19] contain several mathematical results on a chemotaxis system 
and on an aggregation equation with either the fractional Laplacian or a more general 
Levy operator. 

In our recent work [IB] , we have studied the multidimensional version of problem (II. ip - 
(11. 2p and we have answered questions how singularities of the gradient of the aggregation 
kernel K influence on the existence and the nonexistence of global-in-time solutions. In 
this paper, we complete those results in the one- dimensional case by proving the exis- 
tence of global-in-time solutions for every K' £ L X (IR), and by studying their large time 
asymptotics. We show that asymptotic profiles, as t — > oo, of solutions to ( ll.ip - (jl.2p with 
general K' £ L 1 (M) are given either by the fundamental solution of the linear heat equa- 
tion or by self-similar solutions of the viscous Burgers equation. Moreover, under another 
set of assumptions (which are satisfied e.g. by the kernel K(x) = e - ^'), we prove a certain 
concentration property solutions to (ll.ll) - (ll.2p which can be observed on an initial time 
interval. 

Notation. In this work, the usual norm of the Lebesgue space L P (M) is denoted by || • || p 
for any p £ [1, oo] and W k,p (M.) is the corresponding Sobolev space. C^°(M) denotes the 
set of smooth and compactly supported functions. The constants (always independent of 
x and t) will be denoted by the same letter C, even if they may vary from line to line. 
Sometimes, we write, e.g., C = C(a,/3, 7, ...) when we want to emphasize the dependence 
of C on parameters a, 0, 7, .... 

2. Results and comments 

We begin our study of properties of solutions to the initial value problem (ll.ll) - (ll.2l) 
by showing that it has a unique and global-in-time solution for a large class of initial 
conditions and aggregation kernels. The results on the global-in-time existence and reg- 
ularity of solutions from the following theorem are more-or-less standard and we state 
them for the completeness of the exposition. They are also not surprising because it is 
well-known that solutions to the one-dimensional Patlak-Keller-Segel model of chemotaxis 
do not blow up in finite time (see e.g. [251 EE21 II])- On the other hand, the fact that all 
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solutions to problem (ll.lj) — (ll.2p are uniformly bounded in time as the LP- valued functions 
(see Proposition below) seems to be new and is a first step towards the understanding 
of the large time behavior of solutions. 

Theorem 2.1 (Existence of global-in-time solution). Assume that 

(2.1) K' G L\R), 

(2.2) u G L\R) n L q (R) for some 2 < q < oo. 

Suppose also thatu > 0. Then the initial value problem (11.11) - fl 1.2ft has a unique, nonneg- 
ative, global-in-time solution u G C([0, oo), L 1 (IR) fl L 9 (IR)). This solution has the follow- 
ing regularity property u G C 1 ((0, oo), L P {R)) fl C((0, oo), VF 2,P (R))) /or eac/i p G [1, oo]. 
Moreover, the solution conserves the integral ("mass") 

(2.3) M = = / u(x,t) dx = / Uo(x) dx — ||«o||i for all t > 0, 

and for each t > and a//p G (l,oo) we /iave sup t>to < oo. 

In our work [TH], we have proved that nonnegative solutions to (ll.ll) - (ll.2p exist globally- 
in-time, in the case of "mildly singular" kernels satisfying K' G L q (M.) for some q G (1, oo], 
see [TBI, Thm. 2.5]. Theorem 12.11 improves those results in the one dimensional case 
by showing the global-in-time existence of nonnegative solutions to (jl.lj) — (jl.2j) for every 
kernel K' G L 1 (R). This result holds true, in particular, for every "strongly singular" 
kernel pi] satisfying K' G L l {R) \ L q (R) for any q > 1. 

Next, we state conditions on K' under which solutions to (jl.ip - (ll.2j) decay as t — > oo. 

Theorem 2.2 (L p -decay of solutions). Assume that u = u(x,t) is a nonnegative solution 
to problem fll.ll) - fll.2p with K' anduo satisfying (12.11) and (12. 2p . respectively. There exists 
D > and C = C(p, ||i^'||i ||wo||i) > independent oft such that if \\K'\\i\\uo\\i < D, 
then for each p G [1, oo] we have 

(2.4) |K*)IIp<C*~ (1 ~ 1/p)/2 for all t > 0. 

Remark 2.3. In the proof of Theorem I2.2[ we choose D = 1/Cgns, where Cgns is the 
optimal constant in the Gagliardo-Nirenberg-Sobolev inequality (14. 2p . 

We do not know if the decay estimate (12.41) holds true for every nonnegative solution 
which does not necessarily satisfy the condition ll-f^'llill^olli < D. Reasons that such decay 
estimates may fail for certain kernels K' and for initial conditions with large mass can be 
found below in Theorem 12.61 and in the discussion following it. Now, however, we prove 
that estimates (12 ,4p hold true for each solution of problem (ll.ll) - (ll.2p which tends to zero 
as t — > oo without a priori assumed decay rate. 

Theorem 2.4. Let the assumptions of Theorem \2. 1\ hold true. Assume, moreover, that 
there exists po G (1, oo] such that \\u(t)\\ po — > as t — > oo. Then, for each p G [1, oo] there 
is C = C(p, Hif'lji, ||«o||i) > independent oft such that \\u(t)\\ p < Ct~ (1 ~ 1/p)/2 for all 
t > 0. 

The main goal of this work is to derive an asymptotic profile as t — > oo of those solutions 
from Theorem 12. II which satisfy the L p -decay estimates (12.41) . 
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Theorem 2.5 (Self-similar asymptotics). Under the assumptions of Theorem \2.1\ every 
solution u = u(x,t) of problem (11.11) - (11.21) satisfying estimate (12. 4p has a self-similar 
asymptotic profile as t — >■ oo. More precisely, 

i. if J R K'(x) dx = 0, we have 

(2.5) t {1 - 1/p)/2 \\u(t) - MG(t)\\ p ^ as t -4 oo 

for every p G [1, oo], where M = f R Uq{x) dx and Q(x, t) = exp ( — is the 
heat kernel. 

ii. On the other hand, if A = j R K'(x) dx ^ 0, we have 

(2.6) t^- 1/p)/2 \\u{t) -UmA^Wp -> «st^oo 

/or ever?/ p G [1, oo], where Um,a( x ^) = ~Jt^ M > A i^}t' -Q ^ s ^ e s °- ca ^ e d nonlinear 
diffusion wave and is defined as the unique self-similar solution of the initial value 
problem for the viscous Burgers equation 

(2.7) U t = U XX -A (U 2 ) x , for x G E, t > 0, 

(2.8) U(x,0) = M8 , 
where do is the Dirac measure. 

Let us recall properties of solutions to (12.7p - (12.8p which will be useful in the proof of 
Theorem 12.51 It is well-known that the Hopf-Cole transformation allows us to solve this 
initial value problem to obtain the following explicit solution 

(2.9) U M - ^«P(-NV(«» 



Cm,a + Uc , exp(-e/4)^ 

where Cm,a is a constant which is determined uniquely as a function of M and A by the 
condition LUm,a(Vj 1) dr] = M. The important point to note here is that for every MeR 
the function Um,a is a unique solution to equation ( 12 .7p in the space C((0, oo); L 1 (IR)) 
having the properties 



and 



Um,a{ x -> t) dx = M for all t > 



/ ^MA(x,t)(f(x) dx -> M(/?(0) as t -»■ 
Jr 



for all G C£°(R) ( [10l Sec. 4]). Such a solution is called a fundamental solution in the 
linear theory and a source solution in the nonlinear case (cf. [9], [TUJ, |2~T]). 

In the proof of Theorem 12.51 we study the behavior, as A — > oo, of the rescaled family 
of functions 

(2.10) u x {x,t) = Xu{Xx, X 2 t) and K' x {x) = XK'(Xx), for every A > 0, 
which satisfy the initial value problems 

(2.11) d t u x = d 2 x u x - d x (u x (K' x * u x )) , 

(2.12) u x (x, 0) = uq jX (x) = Xuq(Xx). 
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Notice that if u = u(x,t) is a nonnegative solution of problem (jl.ip — (ll.2p obtained in 
Theorem I2.1[ then by (12. 3D and by a simple change of variables, the following identities 

(2.13) IM*)||i=IKI|i and H^Hi = 

hold true for allt > and all A > 0. 

Let us now emphasize that we obtain an asymptotic profile of all solutions to problem 
f)l.ip - f)1.2p which tend to zero as t — > oo. Additional assumptions on the kernel K' and 
on initial conditions (as those in Theorem 12. 2j) seem to be necessary to prove the decay 
of solutions because, in our next theorem, we show a concentration phenomenon for some 
solutions to (jl.ip - (jl.2p . Here, in order to understand our result, one should keep in 
mind that if a solution u = u(x,t) behaves for large t either as the heat kernel or the 
nonlinear diffusion wave, one should expect that Hw^Hoo decreases and the first moment 
J R u(x,t)\x\ dx increases in t. Here, we find a large class of kernels and initial conditions 
such that corresponding solutions have different behavior, at least, on a certain initial 
time interval. 

Theorem 2.6 (Concentration phenomenon). Assume that the kernel K' satisfies 

• K'{x) = K'{\x\) sgnx for all x e R \ {0} ; 

• K'(x) < for all x > 0, 

• there exists 5 > and 7 > such that sup 0<a . <(5 K'(x) < —7. 

Let uq G C£°(R) be nonnegative, nontrivial, and even. For every P > 0, set Mo,p(x) = 
P 3 Uq(Px) and denote by up = up(x,t) the corresponding solution of problem fll.ip - fll.2l) 
with the kernel K' and Wo,p as the initial datum. If P > is sufficiently large, then the 
first moment Ip(t) = f R up(x,t)\x\ dx is a strictly decreasing function of t G [0, T] for 
someT = T(P) > 0. 

Some remarks on Theorem 12.61 are in order. 

Remark 2.7. By the uniqueness, the solution up(x, t) considered in Theorem [22] is an even 
function of x for every t > 0. Hence, it is expected that, for suitable initial conditions, we 
have «p (0,t) = max Ig i«p(i, t). It follows from the proof of Theorem 12 .61 that the quantity 
up(0,t) has to increase on an interval [0, T], however, we do not know if it increases 
monotonically. This phenomenon is in perfect agreement with numerical simulations of 
spikes in the one- dimensional Keller-Segel model, which are reported in |13j . 

Remark 2.8. Assumptions on the kernel K, as those stated in Theorem I2.6[ were imposed 
in [H [16] to show the finite-time blowup of solutions to aggregation equations either in the 
dimension n > 2 or with an anomalous diffusion modeled by the fractional Laplacian. In 
this work, however, solutions are global-in-time by Theorem 12. 11 but they have a tendency 
to form spiky-like structures as those discussed e.g. in [12] . 

Remark 2.9. Notice that M P = J R u 0jP (x) dx = P 2 J R u (x) dx. Hence, the concentration 
phenomenon described in Theorem 12.61 appears only if mass of a solution is sufficiently 
large. This result should be compared with Theorem 12. 2\ where solutions are shown to 
decay for sufficiently small masses. Moreover, in that case, by the inspection of the proof 
of Theorem 12.21 one can show that the L 2 -norm of solutions decays monotonically for all 
t > 0. 
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Remark 2.10. It is shown in Theorem 12 . 61 that the first moment Ip(t) is strictly decreasing 
on a certain initial time interval, however, it cannot converge to zero when t — » oo. Indeed, 
such a decay of Ip(t) cannot be true because of the inequality (see e.g. [3, Remark 2.6]) 

u P (x, t)dx\ < C\\u P (t) \\ p I P (t) l - 1/p , 

the conservation of mass ( 12. 3p . and the boundedness of the L p -norm of solutions to ( ll.ip - 
f ll.2p . shown in Proposition I3.2[ below. 

Remark 2.11. It is not clear for us what is the large time behavior of solutions to problem 
f ll-ip -f TOj) which concentrate initially in the sense described in Theorem 12.61 Our numer- 
ical simulations of solutions to equation (II. ip on a finite interval and with the Neumann 
boundary conditions show their convergence towards nonconstant stationary solutions 
(these results will be published in our subsequent paper). The large time behavior of 
solutions to problem ( ll.ip - (ll.2p on the whole line x G M. seems to be more complicated, 
because one can easily shown that equation (II. ip has no non-zero stationary solutions 
which decay at infinity sufficiently fast. Indeed, the equation w xx — (wK' *w) x = implies 
w x — wK' *w = C, and the constant C has to vanish. Hence, for v(x) = f_ K' * w(y) dy, 
we have (we~ v ) x = and, consequently, w(x) = Ce v<yX \ It is easy to check that w G L X (]R) 
and K' G L 1 (M) implies v G L°°(R), hence, relation w = Ce v can be true for C = and 
w = 0, only. 

Remark 2.12. In [2l| 125] . the authors study the large time behavior of solutions to the 
Cauchy problem for the so-called parabolic-parabolic model of chemotaxis 

(2.14) u t = A«- V • (uVv), v t = Av-v + u, x G IT, t > 0. 

Their results can be summarized as follows. If the solution (u(x,t),v(x,t)) of (I2.14p 
satisfies sup t>0 (||w(t) \\ p + \\v(t)\\ p ) < oo for every p G [l,oo], then u(x,t) decays as a 
solution to the linear heat equation and its large time behavior is described by the heat 
kernel. Moreover, a higher order term of the asymptotic expansion of u is calculated. 
Here, however, because of a technical obstacle, we cannot apply methods from [24"l [25] to 
show a decay of solutions to (ll.ll) - (ll.2p . 

3. Existence of solutions - proof of Theorem 12.11 

The proof of the existence of local-in-time solutions to (ll.ll) - (ll.2p is standard, hence, 
we only sketch that reasoning. 

Step 1. Local-in-time solutions. We construct local-in-time mild solutions of (ll.ll) - (ll.2l) 
which are solutions of the following integral equation 

(3.1) u{t) = Q(.,t)*u - f d x g{-,t- s)* (u{K' *u))(s)ds 

Jo 

with the heat kernel Q(x,t) = (47rt) -1//2 exp ( — |x| 2 /(4t)). In our reasoning, we use the 
following estimates which result immediately from the Young inequality for the convolu- 
tion: 

(3.2) \\G(;t)*f\\ LP <CH(W)||/|| L „ 

(3.3) \\d x g(-,t)*f\\ LP <cr^- l py^\\f\\ L9 
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for every 1 < q < p < oo, each / G L q (R), and C = C(p, q) independent of t, f. Notice 
that C = 1 in inequality (13. 2p for p = q because ||^( , ,t)IU 1 = 1 f° r an ^ > 0. 

Lemma 3.1 (Local existence). Assume that K' G L l (R) and u G L X (1R) n L 9 (IR) for 
some q G [2, oo). Then there exists T = T(||w ||i, ||wo|| g , H-^'lli) > such that the integral 
equation (13.1 p has the unique solution in the space = C([Q, T], L 1 (R))nC([0, T], L q (R)). 
Moreover, this solution satisfies u G C((0, T], L P (IR)) /or all 1 < p < oo. 

Proof. Here, it suffices to follow the reasoning from [MJ Theorem 2.3], where local-in-time 
existence of solutions to the equation (j3.1|) . written as u(t) = <?(•, i) *Mo + -B(u, «)(i) with 
the bilinear form 

(3.4) B(u,v)(t) = - f d x g(;t-s)*(u{K'*v)){s)ds, 

Jo 

are constructed in the space 3^r supplemented with the norm ||u||y T = sup 0<a<T \\u\\i + 
sup 0<t<T \\u\\ g . Notice that <?(-,£) * u G 3^t by (13.21) . To apply ideas from [TBI Theorem 
2.3], one should prove the following estimates of the bilinear form (13. 4p . 

First, for every u,v G using (I3.2p - (I3.3I) and the Young and the Holder inequalities, 
we obtain 

\\B(u,v)(t)\\i < C f{t-s)-V 2 \\u(K'*v){s)\\ x ds 
Jo 

<C I (t-s)-V^u{s)\\M s )\\A K 'h ds i 
Jo 

where 1/g + 1/q' = 1. Since 1 < q' < q for q > 2, by a standard interpolation, we have 
Il v ( s )||g' < + ll'K 5 )!!*?)- Therefore, using the definitions of the norm in 3^r we 

obtain < CT 1 / 2 || J ftr / ||i||w|| 3 ; T ||'y||3; T . 

In a similar way, we prove the following L 9 -estimate 

(3.5) \\B(u,v)(t)\\ q <C [\t-s)- < ^ 1 - 1/q)/2 - 1/2 \\u(K'*v)(s)\\ 1 ds 

Jo 

< ct^iik'ikimi^imi^. 

Summing up these inequalities, we obtain the following estimate of the bilinear form 

||S(tt,«)||^<C7(r 1 / 2 + T 1 /*)||Jf|| 1 ||tt||^||t;||y r . 

Hence, choosing T > such that 4C(T 1 / 2 + T" 1 / 29 ) || AT'H x ( ||m || i + 1Mb) < 1, we obtain 
the solution in y T by [HI Lemma 3.1]. 

Obviously, by an interpolation inequality, the solution of the integral equation (13.11) 
in the space 3V belongs also to C([0, T], L P (R)) for every p G [l,q\- To show that u G 
C((0, T], L P (R)) for all g < p < oo, it is sufficient to note that Q{t) * uq belongs to this 
space by (13.21) . Moreover, for every u,v G 3^t, 

f dx g^ t _ s ) * ( U ( K > * w ))( s ) ds e c((0, T], L p (M)) 
Jo 

for each p G (g, oo] by estimates similar to those in (13. 5p . □ 

Step 2. Regularity of mild solutions. Results on regularity of mild solutions to (II. ip - 
(II. 2p are standard and well-known, see e.g. the monograph by Pazy [27]. In particular, by 
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a bootstrap argument, one can show that any mild solution u G C((0, T], L P (IR)) obtained 
in Lemma O satisfies u G ^((O, T], L P (IR)) n C((0, T], W 2 ' P (M)). 

S^ep 5". Positivity and mass conservation. If the initial condition is nonnegative, the 
same property is shared by the corresponding solution. Obviously, this fact holds true for 
the viscous transport equation u t + u xx + (b(x,t)u) x = 0. Now, it suffices to substitute 
b — K' * u to show that solutions to (11.11) are nonnegative if initial conditions are so. 

Next, integrating equation ( 13. ip with respect to x, using the Fubini theorem, and the 
identities f R G(x, t) dx — 1 and L d x Q(x, t) dx — for all t > 0, we obtain the conservation 
of the integral and the L 1 -norm of nonnegative solutions stated in ( 12. 3p . 

Step 4- Uniform LP -estimate. In order to show that the solution constructed in 
Lemma 13.11 exists for all T > 0, it is sufficient to prove that its L 9 -norm does not blow 
up in finite time. Here, we prove that, in fact, the quantity ||w(i)||p is uniformly bounded 
for large t > for every p G [1, oo). 

Proposition 3.2. Assume that u G C([0, T], L X (R)) n C([0, T],L q (R)) is a nonnegative 
local-in-time solution of problem (ll.ip -f TOj) for some T > 0, with the kernel K' , and the 
initial datum w satisfying (12. ip and ( 12. 2p . respectively. For every p G [1, oo) and for 
every t > 0, there exists C = C{p,to, ||-^'||i, ||wo||<j, H M o||i) independent oft and ofT > 
such that \\u(t)\\ p < C for all t > 0. 

Proof. It follows from Lemma 13.11 that for every to £ (0,7"] we have u(to) G L P (W) for 
all p G [l,oo]. Hence, without loss of generality, we can assume that u G L P {R) and 
t = 0. Moreover, the regularity of solutions discussed in Step 2 allows us to justify our 
calculations below. 

We proceed by induction to show the uniform bound for the L p -norm with p = 2 n , 
n G N PI {0}. The estimates for other p G [1, oo) are a simple consequence of the Holder 
inequality. 

For p — 1, we have = ||«o||i because u is nonnegative and because the integral 

of u is conserved in time, see (I2.3p . 

Now, let e > be small and we fix it at the end of this proof. For p = 2 n with n > 1, 
we multiply equation (II. ip by w p_1 and integrate over M to obtain 



(3.6) 



u p dx 



- / \{u p/2 ) x \ 2 dx + {p- 1) / u p - l u x {K' * u) dx 



pdtJ R p 2 

+ (p - 1) / u p ~ x u x {{K' - K') * uj dx 

J IR. 



where the auxiliary kernel K' G C^°(R) satisfies \\K' — K'\\i < e . 

The second term on the right-hand side of (13. 6p is estimated by the e- Young inequality 
ab < ea 2 + C(e)b 2 and by O as follows 



(p-1) / u p - x u x {K' * u) dx 
Jr 

(3.7) <(p-l)e / u^l^dx + C^IK^ + w)^)!!^ / 



< 4(P 2 1} g / |(« p / a ) x | a tte + C(e)||iiC'||^|KII? / 
P Jr Jr 
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Concerning the third term on the right-hand side of ( 13. 6p . it follows from the identity 
u p ~ 1 u x = (2/p)(u p / 2 ) x u p ' 2 , from the Holder inequality, and from the Young inequality 
that 



(p-1) / u p - x u x {K' -K')*u dx 
2(p-l) 



< 



(3.8) 



< 



< 



P 

2(p-i; 

V 

2(p-i; 

V 



u p, Xu p l 2 \\ P _^\\{K' - K')*u{t)\\ p+ 2 



p+i 

ii(^ /2 yyi^ /2 ik P+2 ) w - K'h\\u(t)\\^ 

V 



2(p+2) > 
V 



because \\K' — K'\\\ < e. Now, notice that the Gagliardo-Nirenberg-Sobolev inequality 
leads to 

IU,P/ 2 II < r\\(,,P/ 2 \ II (P+4)/(3(p+2)) I, p /2||2(p+l)/(3(p+2)) 

1 1 IT' II 2( P +2) < u \\\vr' ) x \\ 2 ||M ' lli 

p 

Moreover, applying the induction hypothesis for p/2 = 2"" 1 , we have ||w p / 2 (t)||i = 
< C with C > independent of t and T > 0. Hence, coming back to in- 



u 



p/2 
p/2 



equality (13. 8p we obtain 



(3.9) 



(p-1) 



vP u x ( (i^' — i^') * it J 



<^IIK /2 ),||f +1)/3p 



with C = C(p, Hi^'Hoo, ||mo||i) independent of t and T. 

Since n(p) = 4(p + l)/3p < 2 for all p > 2, using the elementary inequality s*^ < 
C(s 2 + 1) for all s > and fixed C > independent of s, we deduce from (13.91) that 



(3.10) 



(p-1) 



u p ~ l u x UK' -K')*u) dx 



<Ce\\(u p / 2 ) x \\l + eC. 



Now, by estimates (13. 7p and (13.101) . we obtain from (13.61) the following inequality 



(3.11) 



d 



dt 



- / M ^x<-(l-eL7)||K/ 2 ) x (t)|| 2 + C|| M (t)r + eC, 



where C = C(p, \\K'\\i, ||X'||oo, ||^o||i) denotes various constants independent of e, u, t, 
and T. 

Now, we require e < C~ x in the fist term on the right hand side of (13.111) . Using the 
Nash inequality 

(3-12) \\vh<C\\v x \\y z |M| 2/3 , 

which is valid for all v G L X (IR) such that v x G L 2 (M), and applying the inductive hypoth- 
esis (notice that p/2 = 2" _1 ), we have 

(3-13) |K*)||f = |K /2 (t)|| 2 < C\\ {u p '\t)) x ||f ||^ 2 (t)|| 2/3 < C\\ {u p / 2 (t)) x \\y\ 

where C is independent of u, t, and T. Hence, applying estimate (13.131) in (13. lip we 
obtain the following differential inequality for ||w(t)|| p : 



d_ 
~dt 



W)f v < -C(l - sC)(\\u(t)\\l) 3 + C\\u(t)\\l + eC, 
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where C = C(p, K', ||«o||i) > is independent of e, u, t, and T. 

We leave for the reader the proof that any nonnegative solution of the differential 
inequality /' < —(7(1 — eC)f 3 + Cf + eC is bounded, provided e > is sufficiently small. 
Hence, by the recurrence argument, is bounded for any p = 2 n , neNU {0} and 

this completes the proof of Proposition 13.21 □ 

Remark 3.3. Notice that we do not control the growth in p of the constants C in Propo- 
sition I3.2[ hence we are not able to pass to the limit p — > oo to obtain the global-in-time 
bound of the L°°-norm of the solution. However, in Theorem 12. 2[ we show a decay esti- 
mate of ||ii(£)||oc under additional assumptions on the kernel K' . 

Proof of Theorem \2.1[ Steps 1-4, described above, contain all details of the proof. □ 

4. Optimal L p -decay of solutions 

We are in a position to prove the decay estimates from Theorem 12.21 and we do it 
in two steps. First, we obtain the optimal decay estimate of L 2 -norm using Gagliardo- 
Nirenberg-Sobolev inequalities. Next, estimates of other L p -norms are shown by applying 
the integral formulation (13. ip of the initial value problem ( ll.ip - (11.2p . 

Proof of Theorem \2.2\ Let p — 2. Multiplying equation (11.11) by u and integrating the 
resulting equation over R we have 

(4.1) -— / u 2 dx = — [ \u x \ 2 dx + [ uu x (K' * u) dx 
2 dt J K J R J R 

By the Holder inequality, the Young inequality, the following Gagliardo-Nirenberg-Sobolev 
inequality 

(4.2) \\u(t)\u<c GNS \\uM\l /2 Hml /2 , 

and identity (12.31) . the second term on the right-hand side of (14.11) is estimated as follows 

/ UU x (K'*u)dx< ||^(t)||4/3||(#'*«)(*)l|4 

(43) < IKt)H2K(*)ll2ll^ / lli 

< C , GiV5||Wa ; (*)||2ll^ / ||l||«Q||l- 

Coming back to (14. ip we see that 

^Il«(*)ll2 <-(l-C^5||if||i||«o||i)IM0lla, 
hence, if ||i^'||i||Mo||i < 1/Cgns, we obtain 

(4.4) ^11^(^)111 + ^11^.(^)111 < 0, 

where C = C(\\K'\\i, ||«o||i) > 0. Now, by the Nash inequality (13.121) . since the L 1 -norm 
of the solution is constant in time by (12. 3p . we obtain the differential inequality 

(4.5) |ll^)lll + C||n ||r 4 (||u(t)||l) 3 <0, 

which implies ||u(t)|| 2 < Ct^ 1 ^ for alH > and C > independent of t. 
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Now, we are going to use systematically the above L 2 -decay estimate to show the decay 
of other L p -norms. First, we consider p e [1, oo) and we compute the L p -norm of both 
sides of the integral equation (13.1 j) to obtain 

IK*)||p < \\G(-,t) *Mo|| P + / \\d x Q(-,t-s)*(u(K' *u))(s)ds\\ p ds 

( 4 -6) J ° t 

< C7t~^( 1- *) ||ito||i + C7 / (t-s)"5( 1 -|)-5||^'|| 1 || M ( s )||2 (is 

Jo 

after applying inequalities fl3.2p - fl3.3p . Hence, by the L 2 -decay estimate, we have 

< crK 1 -|)|| Wo || 1 + c'(|| J fr / ||i,||Mo||i) / (t-s)~*(HH s -3d 5 

= Cr^~p) for all t > o, 



o 



because ^(t - s)"^ 1 "^^'^ (is = r^ 1 "?) - p)^ 1- ?) - ^^ dp. 

To deal with the case p = oo, we use the already proved decay estimate for any p > 2. 
Proceeding in a way similar to that in (I4.6p . we obtain 

IK*)||oo<||0(-,*)*«olloo+ / \\d x g(-,t-s)*(u{K'*u)){s)\\ 00 ds 



< Ct~5 ||u ||i + / C(t- s)~p~*\\u(K' *u)(s)\\ p/2 ds 

Jo 

rt 

<c r r3||« ||i+ / c(t-s)"5~3||ir'||iKs)||JcZs 

Jo 

ft 

< cH||« ||i + C(p, H^'lli, IKHO / (*-s)~Hs~( 1_ *)ds 

Jo 

= C(p, \\K'\\ h ||uo||i)H 
for all t > 0. This completes the proof of Theorem I2.2L □ 
In the proof of Theorem 12 A\ we need the following auxiliary result. 

Lemma 4.1. Let v x G L 2 (R). For K' E C~(K), we define K(x) = J^ oQ K'(y)dy. Denote 
A = J R K'(y) dy = lim^+oo K(y). Then 

\\K' *v- AvW^ < (WKW^-oofi) + \\K - A\\ L 2 {0:+oo) j \\v x \\ 2 . 

Proof. First, notice that K e L 2 ((-oo,0]) and K - A e L 2 ([0,+oo)), because K' e 
C£°(R). Hence, the proof is the immediate consequence of the integration by parts and 
of the Schwartz inequality in view of the following identities 

r r+oo 

K' * v(x) - Av(x) = / K(y)v x (x - y) dy - A \ v x (x - y) dy 



□ 



"+0O 

K(y)v x (x - y) dy + I K(y) - A v x (x-y)dy. 
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Proof of Theorem \2.4\ Without loss of generality, we can assume that \\u(t) || 2 — >■ as t — » 
oo, because this is the immediate consequence of the Holder inequality, Proposition 13. 2\ 
and the assumption on the decay of the L Po -norm. 

To show the optimal decay of the L 2 -norm, we use the following equality (cf. (13. 6p ) 



(4.7) 



1 d 



\ u x(t)\\i + / u(t)u x (t)(K' -K') *u(t)dx 
Jr 

I u(t)u x (t)K' *u(t)dx, 



where K e C^°(M) satisfies \\K' — K'\\x < e with e > to be chosen later on. 

It follows from the Holder inequality, the Young inequality, and from the Gagliardo- 
Nirenberg-Sobolev inequality (14. 2 p (see the calculations which lead to (14. 3p ) that 



(4.8) 



u(t)u x (t)(K' - K') * u(t) dx 



< C\\u n 



< Ce||iio||i||u*(*)||l- 



Next, notice that J R u 2 u x dx = for all u G W 1,2 {R), hence, for A = LK'(y) dy, the last 
term on the right-hand side of (14.71) is estimated as follows 



(4.9) 



u(t)u x (t)K' * u(t) dx 



u(t)u x (t) [K' * u(t) - Au(t)) dx 



< \\K' * u(t) - Au^HooHu^) l^llu^,! 112- 

Now, we apply both estimates (14.81) and (14. 9 p in equality (14.71) . Using, moreover, 
Lemma 14.11 to deal with the last term on the right-hand side of (14.91) , we deduce the 
following inequality 

(4.10) 



1 d 
YdV 



u{t)\\ < (-1 + Ce||u ||i + C\\u 



\K\ 



L 2 (-oo,0) 



+ 11^-^41 



L 2 (0,+oo) 



\U X 



Hence, for sufficiently small e > and for sufficiently large To > 0, since \\u 
t — > oo, we obtain the estimate 



2 -> if 



(4.11) 



1 d 
2dt 



\u(t)\\l < -C\\u x (t)\\\ 



for all t >Tq and C > independent of t and u. Now, it remains to repeat the reasoning 
from the proof of Theorem l2.2l for p = 2 (see inequalities (I4.4p - (I4.5I) ) to obtain the required 
decay of the L 2 -norm. To show the decay estimate for other L p -norms, one should copy 
the corresponding arguments from the proof of Theorem 12. 2[ □ 

5. Self-similar large time behavior 

Our goal in this section is to prove Theorem 12.51 Here, we always assume that u = 
u(x,t) is the nonnegative global-in-time solution of the initial value problem (II. !}) - (!! .2p 
with K and Uq satisfying (12 .ip and (12. 2p . respectively. Moreover, we assume that this 
solution satisfies the following decay estimates 

(5.i) \\u(t)\\ p <cr^-7) 
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for each p £ [1, oo], all t > 0, and C independent of t. 

The proof that the large time behavior of the solution u = u(x, t) is described either by 
the fundamental solution of the heat equation or by the self-similar solution of the viscous 
Burgers equation is based on the so-called scaling method which is often use in the study 
of asymptotic properties of solutions to nonlinear evolution equation (see e.g. the review 
article [32] for some applications of this method to the porous media equation). Here, for 
every A > 0, we denote by u x (x,t) = Xu(Xx, X 2 t) the solution of the initial value problem 
fl2.lip - fl2.12p . In the following, we systematically use identities (12.131) as well as the decay 
estimate (15. ip . 

Now, we prove a series of technical lemmas which usually should be obtained to apply 
the scaling method. 

Lemma 5.1. For each p £ [1, oo] there exists C = C(\\K'\\i, ||wo||i) > ; independent of 
t and of X, such that 

(5.2) KWIIp^crK 1 -*) 

for all t > and all X > 0. 

Proof. By the change of variables and estimate (15. ip we obtain 

\\ Ux (t)\\ p = x 1 - 1 ? ik, x 2 t)\\ p < cx 1 - 1 ? (A 2 tp (l ~^ = cHO-l). 

□ 

Lemma 5.2. For each p £ [1, oo) there exists C = C(p, \\K'\\i, ||wo||i) > 0, independent 
oft and of X, such that \\d x u\(t)\\ p < Ct~^~^'^ for allt > and all X > 0. 

Proof. Here, we use the the following counterpart of the integral equation ( 13. ip 

d x u x (t + i) = d x g(t)*u x (i) 

- j d x Q{t - s) * {{u x ) x (K' x * u x ) + it a (K' x * (u x ) x ) ) (s + 1) ds 
for all t > 0. Hence, computing the L p -norm and using (13. 2p . (13. 3p . (I2.13P we obtain 

II^c«a(*h- i)||p < ct-3( 1 -F)-i|j« A (i)|j 1 

+ J C(t - s )-5 (H (u x ) x (K' x * u x ) (s + 1) \\ p + \\u x {K' x * (u x ) x ) (s + 1) || p ) ds 

+ / (t-s)-^\\d x u x (s + l)\\ p \\u x (s + 1)^ ds. 

Jo 

Next, we use inequality (15.21) with p = oo: 

\\u x (s + l)!^ < Cdl^'Hx, ||« ||i)(s + I)"* < Cdl^'Ha, ||« ||i) 
for all s > 0, consequently, 

Wd^t + vWp^cr^-^Wuoh + c [ (t- s )-^\\d x u x (s + i)\\ p ds. 

Jo 
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Applying the singular Gronwall lemma (see e.g. [Hj Ch. 7] ), we conclude that 

+ < ^feM^'lli, IKIli), 

for all t > 0, where the right-hand side of this inequality is independent of A. In particular, 
for t — 1, there exists C = C(\\K'\\i, ||wo||i) independent of A such that \\d x u\{2) \\ p < C for 

all A > 0. Next, using the definition of u\, we obtain ||(9 :e -ua(2)||p = X 2 ~p\\d x u(-,2X 2 )\\ P . 

Hence, after substituting A = yt/2, we arrive at \\d x u(-,t)\\ p < (7r*( 1_ *H for all 
t > 0. □ 

The proof of our next lemma relies on a form of Aubin-Simon's compactness result that 
we recall below. 

Theorem 5.3 ([28, Theorem 5]). Let X , B and Y be Banach spaces satisfying X C B C 
Y with compact embedding X C B. Assume, for 1 < p < oo and T > 0, that 

• F is bounded in L p (0,T; X) , 

• {d t f : / eF} is bounded in 17(0, T;K). 

Then F is relatively compact in L p (0,T; B) (and in C(0,T; B) if p = oo). 

Lemma 5.4 (Compactness in Lj oc (M.) ). For every < t\ < ti < oo and every R > 0, 
the set {«a}a>o — C([^i> ^2], L 1 ^— R, R])) is relatively compact. 

Proof. We apply Theorem 15.31 with p = 00, F = {u\}\ >0 , and 

X = W x,1 {[—R, R)), B = L\[-R, R)), Y — W~ 1,l {[—R, R}), 

where R > is fixed and arbitrary, and Y is the dual space of Wo' ( [-R, R] ) . Obviously, 
the embedding X C B is compact by the Rellich-Kondrashov theorem. 

By Lemmas O and El with p = 1, the sets {w A } A>0 C L°° ([t 1} t 2 ], ^([-R, R})) and 
{dsudx^CL^ifahlL^-R, R})) are bounded. 

To check the second condition of Aubin-Simon's compactness criterion, it is suffices 
to show that there is a positive constant C which independent of A > such that 
sup t6 r tl]f2 i Hc^maHy < C. Let us show this estimate by a duality argument. For every 
(f) E ((—R, R)) and t e [t!,t 2 ] we have 

d t u\(x, t)4>(x) dx = — I d x u\d x (f)dx + / (d x 4>)u\(K' x * u\) dx 

Jr Jr 

< ||^||oo(||UAWIIl+ll^lllH«A(t)||l) 

< H0xllooC(ti,t 2 , ll^'lli, IMIO 

by virtue of Lemma [5.1[ Hence, Lemma [5.41 is proved. □ 

Lemma 5.5 (Compactness in L 1 (IR) ). For every < t\ < £ 2 < 00, the set {ma}a>o Q 
C([ti, t 2 ], L 1 (IR)) is relatively compact. 

Proof. Let if) G C°°(1R) be nonnegative and satisfy ip(x) =0 for \x\ < 1 and ip(x) = 1 for 
\x\ > 2. Put i/jr(x) = ip(x/R) for every R > 0. Since u is nonnegative, in view of Lemma 
15. 4\ using a standard diagonal argument, it suffices to show that 

(5.3) sup ||ma(^)V'h||i — > as R — > 00, uniformly in A > 1. 

te[ti,t 2 ] 



ONE-DIMENSIONAL MODEL OF CHEMOTAXIS 



15 



Multiplying the both sides of equation ( 12.1 ip by ip R and integrating over R and from 
to t, we obtain 

ux(x,t)ip R (x) dx — / u\{x, 0)ijj R (x) dx 





d xx ip R (x)u\(x, s) dxds + d x tp R (x) (u\(x, s) (K' x * u x ) (x, s)) dxds. 



Since d xx ip(x) = ip" (x / R) / R 2 and d x ip R (x) = \x / R) / R, we have 



u\(x, t)ijj R (x) dx < u\fl{x)il) R {x) dx + 



R 2 



||n A (s) ||i ds 



+ 



(5.4) 



R 



< / ux,o(x)ip R (x) dx 



WKh / IKOOIUM^M* 



'lloolKlllll^'lh 



+ -^|kMIoo||uo||i + 

because ||wa(s)||i = ||uo||i- 

Now, notice that, by the change of variables, we have 



R 



||ma(s)||oo ds 



u\fl(x) dx 



\x\>R 



uo{x) dx < 



\x\>XR 



uq(x) dx 



\x\>R 



for all A > 1. Moreover, it follows from Lemma 15.11 that ||wa(s)||oo < Cs l l 2 with C 
independent of A. Hence, we see for all t G [ti,t 2 ] that 



u\{x,t)ip R (x) dx < 



\x\>R 



t 2 t 1/2 ' 

u {x) dx + C 2 | + 



where C 2 = C^GlVksllooj Halloo, ll-^'lli) |l w o||i) is independent of A > 1. This proves our 
claim (15 .3p because uq G L 1 (M). □ 

Lemma 5.6 (Initial condition). For every test function <p G C£°(R), there exists C = 
C(4>, \\K'\\i, || Mo ||i) independent of X such that 



(5.5) 



u\(x,t)(j)(x) dx — / uo,x(x)(j)(x) dx 



KC^ + t 1 ' 2 ) 



Proof. Following the calculations which lead to estimates (15. 4 p with ip R replaced by G 
C°°(R) we obtain 



u\{x,t)(j)(x) dx — I uo,x(x)(j)(x) dx 



< 



\\ux(s)\\ids + 



l\K'MMs)\\i\\Ms)\\oods 



< ||0xx||oo||«o||l* + C1l<Ar||oc||^'||l |l"<i|l I 



<C(t + t 1/2 ), 
where C > is independent of A. 



□ 
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Now, we are in a position to prove our main result on the large time behavior of solutions 
to problem (O]) - (T01) . 

Proof of Theorem \2.5[ By Lemma I5"3l for every < t\ < t 2 < oo, the family {wa}a>o 
is relatively compact in C([ti, t 2 ], L 1 ^)) for any < ti < t 2 < oo. Consequently, there 
exists a subsequence of {wa}a>o ( n °t relabeled) and a function u G C((0, oo), L 1 (M)) such 
that 

(5.6) u x u in C([t u t 2 ], ^ W) as A ^ oo. 
Passing to a subsequence, we can assume that 

(5.7) u\(x, t) — )• u(x, t) as A — > oo 

almost everywhere in (x, t) £ 1 x (0, oo). 

Now, multiplying equation (12. lip by a test function e C c °°(lx (0, oo)) and integrating 
the resulting equation over R x (0, oo), we obtain the identity 



(5.8) ~ u\4> t dsdx = I I u\(j) xx dsdx + / / U\{K' X * u\)4> x dsdx. 






Recall that K' x (x) = \K\\x) and A = f R K'(y) dy = J R K' x (y) dy for all A > 0, hence, 
by the well-known property of an approximation of the identity, we have K' x *u(t) — > Au(t) 
in L X (IR) as A — > oo. Consequently, by the Young inequality, (12 . 1 3j) . and (15. 6p . we have 

\\K' X * u x (t) - Auit)^ < WK'WxWuxit) - tZ(t)||i + ||^ * u(t) - Auit)^ -)• 

as A — > oo for every t > 0. Therefore, passing to the limit A — > oo in equality (15. 8p and 
using the properties of the sequence {«a}a>o stated in (I5.6P and (15. 7p . we obtain that 
u(x,t) is a weak solution of the equation 

u t = u xx - A{u 2 ) x 

with A = J R K'(x) dx. Now, notice that, by the change of variables and the dominated 
convergence theorem, we obtain J R uo t \(x)(p(x) dx = J R uo(x)(f)(x/ A) dx — > M0(O) as A — > 
oo. Hence, it follows from Lemma 15.61 that u(x, 0) = M5q in the sense of bounded 
measures. Thus, u is a weak solution of the initial value problem 



(5.9) u t = u xx - A(u Jx , 

(5.10) u(x,0) = M5 . 

Since problem (l5.9p -f l5TTU|) has a unique solution (see e.g. [TUJ Sec. 4]), the whole family 
{^a}a>o converge to u in C((0, oo), L 1 (M)). 

Obviously, if A = 0, this limit function is the multiple of Gauss- Weierstrass kernel 

1 Id 2 

(5.11) u(x,t) = MG(x,t) = M-^=exp(- LL). 

For A 7^ 0, we obtain the self-similar solution u = Um,a of the viscous Burgers equation, 
given by the explicit formula (12. 9p . 
Hence, by (I5.6p . we have 

lim ||« A (l)-«(l)|| a = 

A— >oo 
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and, after setting A = y/t and using the self-similar form of u(x,t) = t~ l l 2 u(xt~ l l 2 , 1), we 
obtain 

(5.12) lim \\u{t) -u(t)||i = 0. 

t— >oo 

The convergence of u(-, t) towards the self-similar profile in the L p -norms for p e (1, oo) 
is the immediate consequence of the Holder inequality, the decay estimate 05. ip with 
p = oo, and (I5.12p . Indeed, we have 

(5.13) \\u(t) - «(t)|| p < (||«(t)|U + II^^IU) 1 ^^^) - u(t)\\l /p = o^ 1 - 1 ^ 2 ) 
as t — > oo. 

To complete the proof of Theorem 12. 5[ it remains to show the convergence in the L°°- 
norm. Here, however, it suffices notice the decay estimate ||ua;(t)||2 < Ct~ 3 / 4 , provided 
by Lemma I5T21 with A = 1, and the identity ||u x (t)||2 = ^ 3 ^ 4 ||^x(l)||2 resulting from the 
explicit formulas (15. lip and (12.91) . Hence, by the Gagliardo-Nirenberg-Sobolev inequality 
and by (15.131) with p = 2, we obtain 

\\u(t) - w^lU < C{\\u x (t)\\ 2 + \\u x (t)\\ 2 ) 1/2 \\u(t) - u(t)\\l /2 = o(t- 1 ' 2 ) 
as t — > oo. □ 

6. Concentration phenomenon 

Proof of Theorem \2.b\ First, notice that equation equation (II. ip is invariant under the 
transformation x \- > —x because, by the assumptions, the kernel K' is an odd function. 
Hence, by the uniqueness of solutions to problem (ll.ip - (ll.2l) . the solution up corresponding 
to the even and nonnegative initial datum uo,p satisfies up(x,t) = up(—x,t) > for all 
x e R and t > 0. 

We study the evolution of the first moment Ip(t) = L Up(x, t)\x\ dx following ideas 
from the recent papers [H [16]. First, notice that integrating by parts we have 

p poo i*oo 

/ \x\d 2 up{x, t) dx = 2 / xd 2 up(x, t) dx = — 2 / d x up(x, t) dx = 2up(0, t), 
Jr Jo Jo 

because d x up(0, t) = for alH > in the case of the even function up(-,t). Next, by the 

assumptions, we have K'{x) = ■AK'(\x\) where ^ = sign x. Hence, using equation (11.11) . 

we obtain 

d f 

—I P (t)= I d t u P (x,t)\x\dx 

OX J TO 



(6.1) 



(d 2 up(x, t) — d x (up(x, t)K' * up(x, t))) \x\ dx 

f f xx — y 

= 2up(0,t)+ / / up(x,t)up(y,t)- — 1| \K'{\x — y\) dxdy 

Jr Jr fI F — HI 

1 f f fx y \ x — y 

= 2u P (0, t) + - / u P (x, t)u P (y, t) I | — | — -j — | I | ;K'(\x - y\) dxdy 

2 Jr Jr VfI \y\J \x-y\ 

by the symmetrization of the double integral on the right-hand side. Now, we use the 
elementary identity 

x y\ x — y f\x\ + \y\\ ( x y 



x\ \y\j \x — y\ \ \x — y\ j \ \x\ \y 
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and the inequalities 

\x\ + \v\ , xy 

1 < V and 1 - -- ■ f- > 

\x — y\ \x\ \y\ 

which are valid for all x,y G R \ {0}. Moreover, using the properties of K' stated in 
Theorem 12 .6[ we deduce from (16. ip that 

j t I P {t) < 2u P (0, t) + l - J J u P (x, t)u P (y, t) (l - • Hfflx - y|) dxdy 



M«5/2 
|w|«5/2 



< 2w P (0,t) - I yy Up(x,t)u P (y,t) ^1 - • ) f/.rJ/y 



(6.2) i^i<«/2 

l»l<«/2 



< 2u P (0,t) 
+ 7 





1 2/1x5/2 





a; 






I 1 


F 


a; 




;r 


\y\j 



Note now that 



fx fx 
/ Mp(x,t)-j — | dx = and / Mp(i, t) — - = 

JR fI J\x\>5/2 Fl 



|x|>.5/2 

because up(x,t) = up(—x,t). Hence, denoting Mp = J R up(x,t) dx it follows from (I6.2p 
that 

-I P {t) < 2u P {0, t) - ]-M 2 P + 7 M P / u P (y, t) dy 

(6.3) 



\y\>&/2 



7 ,.2 , 2 7 



< 2« P (0,t) - + -±M P Ip{t). 
2 o 

Now, we use the dependence of the initial datum on P. Putting 

M= u (x)dx, 1(0) = / u (x)\x\dx, 



and changing the variables we have 

2« p(0) - + ^M P /p(0) = 2m (0)P 3 - ^M 2 P 4 + ^M/(0)P 3 < 
2d 2d 

if P > is sufficiently large. Hence, for such P, the right hand side of inequality ( 16. 3p is 
negative for t = . Consequently by the continuity of the functions up and Ip, the right 
hand side of inequality (16. 3p is negative for every t G [0,T] with some T = T(P) > 0. 
This completes the proof of Theorem 12.61 □ 

References 

[1] A.L. Bertozzi, J. A. Carrillo, and T. Laurent, Blowup in multidimensional aggregation equations with 

mildly singular interaction kernels, Nonlinearity 22 (2009), 683-710. 
[2] A.L. Bertozzi and T. Laurent, Finite-time blow-up of solutions of an aggregation equation in R", 

Comm. Math. Phys. 274 (2007), 717-735. 
[3] P. Biler, G. Karch, Blowup of solutions to generalized Keller-Segel model, J. Evol. Equ. (2010), 

published online 6 November 2009. 



ONE-DIMENSIONAL MODEL OF CHEMOTAXIS 



19 



P. Biler, G. Karch and Ph. Laurengot, Blowup of solutions to a diffusive aggregation model, Nonlin- 
earity 22 (2009), 1559-1568. 

A. Blanchet, J. A. Carrillo, Ph. Laurengot, Critical mass for a Patlak-Keller-Segel model with degen- 
erate diffusion in higher dimensions, Calc. Var. Partial Differential Equations 35 (2009), 133-168. 
A. Blanchet, J. Dolbeault, B. Perthame, Two dimensional Keller-Segel model: Optimal critical mass 
and qualitative properties of the solutions, Electron. J. Diff. Eqns. 2006, 44, 1-33. 
M. Bodnar, J.J.L. Velazquez, Derivation of macroscopic equations for individual cell-based models: 
A formal approach, Math. Meth. Appl. Sci. 28 (2005), 1757-1779. 

M. Bodnar, J.J.L. Velazquez, An integro- differential equation arising as a limit of individual cell- 
based models, J. Differential Equations 222 (2006), 341-380 

I.L. Chern and T.P. Liu, Convergence to diffusion waves of solutions for viscous conservation laws, 
Comm. Math. Phys. 110 (1987), 1103-1133. 

M. Escobedo, J.L. Vazquez, E. Zuazua, Asymptotic behaviour and source-type solutions for a 
diffusion-convection equation, Arch. Rational Mcch. Anal. 124 (1993), 43-65. 

C. Escudero, The fractional Keller-Segel model, Nonlinearity 19 (2006), 2909-2918. 
T. Hillen, A classification of spikes and plateaus, SIAM Review 49 (2007), 35-51. 

T. Hillen and A. Potapov, The one- dimensional chemotaxis model: global existence and asymptotic 
profile, Math. Methods Appl. Sci. 27 (2004), 1783-1801. 

D. Henry, Geometric Theory of Semilinear Parabolic Equations, in: Lecture Note in Mathematics, 
Vol. 840, Springer, Berlin, 1981. 

T. Ikeda, T. Nagai, Stability of localized stationary solutions, Japan J. Appl. Math. 4 (1987), 73-97. 

G. Karch and K. Suzuki, Blow-up versus global existence of solutions to aggregation equation with 
diffusion, (2009), l-14. iarXiv:1004.4021l 

H. Kozono and Y. Sugiyama, Local existence and finite time blow-up of solutions in the 2-D Keller- 
Segel system, J. Evol. Equ. 8 (2008), 353-378. 

T. Laurent, Local and global existence for an aggregation equation, Commun. Partial Diff. Eqns 32 
(2007), 1941-1964. 

D. Li and J. Rodrigo, Finite-time singularities of an aggregation equation in M™ with fractional 
dissipation, Comm. Math. Phys. 287 (2009), 687-703. 

D. Li and X. Zhang, On a nonlocal aggregation model with nonlinear diffusion, Discrete Contin. Dyn. 
Syst. 27 (2010), 301-323. 

T.P. Liu and M. Pierre, Source solutions and asymptotic behavior in conservation laws, J. Differential 
Equations 51 (1984), 419-441. 

J. D. Murray, Mathematical Biology. II, Spatial Models and Biomedical Applications, vol. 18 of 
Interdisciplinary Applied Mathematics, Springer- Verlag, New York, third ed., 2003. 
T. Nagai, Behavior of solutions to a parabolic- elliptic system modelling chemotaxis, J. Korean Math. 
Soc. 37 (2000), 721-732. 

T. Nagai, R. Syukuinn, M. Umesako, Decay properties and asymptotic profiles of bounded solutions 
to a parabolic system of chemotaxis in R N , Funkcial. Ekvac. 46 (2003), 383-407. 
T. Nagai, T. Yamada, Large time behavior of bounded solutions to a parabolic system of chemotaxis 
in the whole space, J. Math. Anal. Appl. 336 (2007), 704-726. 

K. Osaki, A. Yagi, Finite dimensional attractor for one- dimensional Keller-Segel equations, Funkcial. 
Ekvac, 44 (2001), 441-469. 

A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, 
Springer- Verlag, New York, Berlin, Heidelberg, Tokyo, 1983. 

J. Simon, Compact sets in the space L p (0,T; B), Ann. Mat. Pura Appl. 146 (1987), 65-96. 

A. Stevens, A stochastic cellular automaton modeling gliding and aggregation of myxobacteria, SIAM 

Journal on Applied Mathematics 61 (2000), 172-182. 

Y. Sugiyama, On e-regularity theorem and asymptotic behaviors of solutions for Keller-Segel systems, 
SIAM J. Math. Anal. 41 (2009), 1664-1692. 

C. M. Topaz, A. L. Bertozzi and M. A. Lewis, A nonlocal continuum model for biological aggregation, 
Bull. Math. Biol., 68 (2006), 1601-1623. 

J.L. Vazquez, Asymptotic behaviour for the porous medium equation posed in the whole space. Ded- 
icated to Philippe Benilan, J. Evol. Equ. 3 (2003), 67-118. 



20 



GRZEGORZ KARCH AND KANAKO SUZUKI 



Instytut Matematyczny, Uniwersytet Wroclawski, pl. Grunwaldzki 2/4, 50-384 Wroc- 
law, POLAND 

E-mail address: grzegorz.karch@math.uni.wroc.pl 
URL: http : //www .math. uni .wroc .pl/~karch 

Institute for International Advanced Interdisciplinary Research, Tohoku University, 
6-3 Aramaki-aza-Aoba, Aoba-ku, Sendai 980-8578, JAPAN 
E-mail address: kasuzu-is@rn.tains.tohoku.ac.jp 



